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Abstract 

Thin ferromagnetic elements in the form of rectangular prisms are theoretically 
investigated in order to study the transition from single-domain to two-domain state, 
with changing the in-plane aspect ratio p. We address two main questions: first, how 
general is the transition; second, how the critical value p c depends on the physical 
parameters. We use two complementary methods: discrete-lattice calculations and 
a micromagnetic continuum approach. Ultrathin films do not appear to split in two 
domains. Instead, thicker films may undergo the above transition. We have used the 
continuum approach to analyze recent Magnetic Force Microscopy observations in 
30nm-thick patterned Permalloy elements, finding a good agreement for p c . 
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1 Introduction 

In recent years, arrays of patterned ferromagnetic dots have received consid- 
erable interest owing to their possible applications in high-density magnetic 
data storage [1] and spin-electronic [2] devices, as well as for realizing logic 
functionality [3]. High- resolution electron beam lithographic techniques [4] are 
commonly used to fabricate the samples, in such a way that all the particles in 
the array are virtually identical to each other. As a consequence, the measured 
properties of the array can be interpreted as the individual properties of a sin- 
gle dot, provided that the dots are sufficiently far spaced to neglect the interdot 
magnetostatic interaction. In this way, using high-sensitivity magneto-optical 
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magnetometry techniques, the variation of the properties with the shape in 
magnetic nanoelements could be experimentally investigated for dots as thin 
as 3 nm [5]. Single dot properties can also be investigated by magnetic force 
microscopy (MFM) techniques [6], which allow to focus on a single element. 

In this paper, we will consider in-plane magnetized dots [1]. Within the plane, 
the preferred direction of the magnetization is determined by a balance be- 
tween magnetic dipole-dipole interaction and magnetocrystalline anisotropy [7]. 
In Permalloy elongated nanodots, the direction of the uniform magnetization 
is determined by the shape-induced magnetic anisotropy, since the crystalline 
anisotropy of the material is negligible. 

In a recent paper [8] we studied rectangular monolayers of planar spins lo- 
cated on the sites of a two-dimensional triangular lattice and interacting via 
magnetic dipole-dipole interaction only. For this lattice geometry, the infinite 
monolayer (ML) is believed [9] to have a ferromagnetic (FM) ground state, 
which is degenerate with the orientation. The infinite square lattice is instead 
believed to have a microvortex ground state, which is degenerate with a local 
orientation angle. For specific values of this angle, we get a state of FM lines 
antiferromagnetically coupled. Therefore, the difference between the two types 
of lattice is often understood in terms of the coupling between FM lines of 
spins, which is ferromagnetic in the triangular lattice and antiferromagnetic 
in the square lattice. 

In our paper [8], rectangular monolayers were studied with varying their size 
and aspect ratio p, defined as the ratio between the smaller and the larger side 
of the rectangle. 1 Two main results came out. First, passing from rectangular, 
elongated finite systems to square-shaped ones, we obtained clear evidence for 
a transition from a FM to a two-domain and eventually to a macrovortex 
state configuration. Second, in the thermodynamic limit, we found that the 
macrovortex seems to be the lowest energy state whatever is the aspect ratio, 
but for elongated samples the size must attain unphysically large values to 
display such a state. 

In this paper, we want to study the effect of a direct exchange interaction 
between spins, taking into account the effect of a variable thickness as well. 
These generalizations make our result relevant for FM prisms which are cur- 
rently available. A recent paper [11] by Jubert and Allenspach moved along 
similar lines, with the authors studying the transition from single-domain 
to macrovortex configuration for a circular ferromagnetic disk. Here, we are 
rather interested in the problem of the transition from single-domain to two- 

1 Such a definition was adopted here for better convenience as regards the com- 
parison with the experimental data in Ref. [3]. Note that it differs both from the 
one adopted by Aharoni in Ref. [10] and from the one we used in Ref. [8] (r = l/p 
rather than p) . 
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domain configuration. For this reason, we study a thin rectangular ferromag- 
netic prism. 

We want to address two main questions: 

i) How general is the transition? We have already shown [8] that this transition 
occurs for a purely dipolar system in the ML limit. Introducing the exchange 
coupling, A, makes the domain wall more energetically expensive and change 
the width of the domain wall: the latter is atomically sharp for a purely dipolar 
system, while it spreads over many lattice constants when A ^ 0. Here, we 
want to ascertain whether these changements are able to cancel the transition 
and whether lattice structure effects, which are so important for a purely 
dipolar system, are maintained. 

ii) If a transition occurs, how does the critical value, p c , of the in-plane aspect 
ratio, p, of the rectangular element depend on the physical parameters of the 
system? 

We will use two complementary methods for the evaluation of the magne- 
tostatic energy difference, AEm, between a two-domain state and a single- 
domain state: discrete-lattice calculations, which are valid for ultrathin films, 
and a micromagnetic continuum approach, applicable to thicker films. 

In Section II discrete-lattice calculations are performed for a rectangular mono- 
layer, both for spins located on a triangular and on a square lattice. Our main 
result is that the transition from a single-domain state to a two-domain state 
is suppressed, in the ultrathin limit, because of the exchange-induced broad- 
ening of the domain wall, which modifies the energetic balance of the dipolar 
energy. This result does not depend on the type of lattice, because spin con- 
figurations varying on scales much larger than the lattice constant depend 
weakly on lattice geometry. 

In the case of thicker samples, a micromagnetic approach is appropriate. In 
Section III we calculate the magnetostatic energy difference, AEm, between 
the two-domain and the single-domain state, following a method devised a 
few years ago by Aharoni [10]. In the limit of domain-wall width much smaller 
than the lateral dimensions of the ferromagnetic dot, analytical expressions 
can be obtained both for the magnetostatic energy gain, AEm, due to surface 
charges and for the energy cost (per unit wall area), 77V, of a one dimensional 
Neel wall [12,13]. For a thin rectangular prism, 7^ includes contributions from 
the exchange interaction, the uniaxial magnetocrystalline anisotropy, and the 
magnetostatic energy due to volume charges in the Neel wall. One can thus 
estimate in a simple way the critical value, p c , of the in-plane aspect ratio 
separating the single-domain phase from the two-domain phase. 

An important remark is in order here. Introducing the exchange interaction, 
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A, has a twofold effect: i) the domain wall width is increased to a value L& w ss 



yjA/Ki, where K\ is the uniaxial anisotropy constant, and ii) the domain wall 
energy density (i.e., the energy per unit wall area) is increased by a quantity 
of order \J AK\. A two-domain state can be energetically favorable if the gain 
in dipolar energy density coming from the two opposite domains, of order M^b 
(where M s is the saturation magnetization and 2b is the thickness), prevails on 
the full domain wall energy density, which has a dipolar contribution plus the 
just mentioned term of order \[AK[. A necessary condition for the transition 
to occur is that the two-domain configuration is energetically favorable when 
A = 0, because the exchange (and the anisotropy) only contribute to the 
cost of the energy balance, not to the gain. If the transition does occur when 
A = 0, one can wonder whether the transition is maintained when A ^ 0. 
One is thus led to conclude that, for comparable values of the anisotropy and 
dipolar energy densities (Ki ~ M%), the transition disappears because the 
domain wall energy density (of order y/AIC[) surely dominates on the gain 
(of order M^d). Such a no-domain rule for in-plane configurations in ultrathin 
films had already been found in Refs. [7,14]. 

However, there are systems where K\ is so small that the term \[AK\ is 
comparable to, or even smaller than, M^b: this is just the case of Permalloy. 
In fact, our theoretical predictions about the occurrence of a transition from 
a single-domain to a two-domain state in rectangular Permalloy dots when 
increasing the in-plane aspect ratio p are confirmed by recent MFM data 
[3] of elliptic patterned Permalloy elements with thickness 30 nm, lateral size 
3 jum and p varying between the values p — 0.13 (elongated dot, single-domain 
state) and p < 1 (almost circular dot, two-domain state). 

Finally, in Section 4 we compare the discrete and the continuum approaches. 
There are some indications that the transition from single-domain to two- 
domain state may appear when the thickness of the film is comparable to the 
width of the domain wall. Final considerations are reported in Section 5. 



2 Discrete-lattice calculation of the magnetic dipolar energy of a 
one-monolayer-thick rectangular element 

Let us start considering a rectangular element which is one monolayer thick. 
We want to compare the energy of the in-plane ferromagnetic (FM) configu- 
ration, where all spins (because of shape anisotropy) are parallel to the longer 
side, to the two-domain configuration (see Fig. 1), where a single Neel do- 
main wall, parallel to the longer side, is inserted in the middle of the stripe. 
We have considered both a triangular and a square lattice, with the z axis 
parallel to a row of atoms. However, for ease of notations, we will refer to 
the square lattice. Different rows are labelled with index m as follows: for 
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Fig. 1. The geometry of the rectangular element in the two-domain configuration. 
The bold arrows denote the magnetization directions and 2Ld w is the domain wall 
width. 

—L x < m < — 1 (region '1'), the spins are parallel to the +z direction, S = z; 
when m varies between m = and m = 2Ld w (region '0'), the spins rotate 
inside the domain wall, according to the relations S x (m) = cos[m7r / '{2L^J)\ 
and S y (m) = sm[mn / (2L dw )]); finally, for 2L dw + 1 < m < 2L dw + L x (region 

— * 

'-1'), S = —z. If the distance between nearest neighboring spins is taken as 
length unit, the size of the stripe in the z direction is 2L Z and lines parallel 
to z are distant 1. 

The full dipolar energy of a given configuration can be written as 



(ij) Ti i 



S{ ' Sj 3 



(Si ■ Tij)(Sj ■ fij) 



where Q = fi 2 /al = M^a^ (/i is the magnetic moment per spin, M s is the 
magnetic moment per unit volume, and ao is the lattice constant). 

Interacting spins (ij) can be grouped according to the regions, '1', '0' and '-1', 
they belong to, so that: 

Edi P = Ei t i + £o,o + £-1,-1 + -£1,0 + £-1,0 + Ei -1. (2) 



The first three terms are the self-energies of the three regions; Ei fl and E_ 1>0 
are the interaction energies between each ferromagnetic region and the do- 
main wall region; finally, £1,-1 is the energy between the two ferromagnetic 
regions. We are interested in the energy difference between the two-domain 
state and the FM (single-domain) state. Since it will be useful to consider ener- 
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gies per unit length in the z direction, we define AE dip = (E dip — E di ¥) / (2L Z ), 
where E dip and are the energies of the two-domain and of the FM 

state, respectively. When the symbols E and E FM are accompanied by the 
subscripts '1', '0', '-1', the energies refer to self /interaction dipolar energies 
of specific regions. Since regions '1' and '-1' keep themselves ferromagnetic, 
A#i,i = A£_i_i = 0, and 

AE dip = AE , + AE lfi + AE_ lfi + Af?i _i. (3) 



As for AE^ and AE_i t0 (which are equal, for symmetry reasons), they can 
be neglected [15,16], because the interaction between a FM domain and the 
domain wall is averaged to zero in the limit L x ^> L dw . Finally, we have 

AE ~ A£ ,o + AE 1 ! _i, (4) 



where the first term is the self-energy of the domain wall and the second term 
is the interaction energy between the two FM regions. In Fig. 2 (main) we 
show that AEofl goes to a constant with increasing L z , so that we can define 
a dipolar domain wall energy per unit length. As for the L^-dependence of 
AE 0t0 , it is shown in the inset of the same figure. For a large range of values 
of L dw , 1 <C L dw <C L z , AE 0j q can be approximately taken as constant. In 
conclusion, we can assume that AE 0t0 /Q is a constant do, with do — 4. 

The quantity A£i ; _i can be rewritten as follows, 

A-Ei ,-i = (iVi - E™)/(2L Z ) = -2E™/(2L Z ), (5) 



because in the two-domain state, the spins in the '-1' region are just reversed 
with respect to the FM state. Since the two regions are 2L dw far away, we can 
replace the discrete summation by an integral. We get the result: 



2sju? z + 2L 2 dw - A^ALl + (L x + 2L dw f + 2^L\ + (2L X + 2L dw f 
— ^{L x + 2L dw ) \n(L x + 2L dw ) + 2L dw In 2L dw 



- (QL X + 8L dw ) ln[^4L2 + (L x + 2L dw f — L x — 2L dw ] 

- 2L X In [y/AL* + (L x + 2L dw f + L x + 2L dw ] 

+ 2(2L X + 2L dw ) \n(2L x + 2L dw ) + 2(L X + 2L dw ) H^AL* + 2L\ W - 2L dw ] 
+ + 2L dw \ 

+ 2(2L X + 2L dw ) H^ALl + (2L X + 2L dw f - 2L X - 2L dw }. (6) 
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In the relevant limit Ld w <C L Z ,L X , we get 



AEi_i/fi ~ -4^ 1 + pln4 - \/4 + p 2 + \/l +p 2 - pin 



( 



a/1 + p 2 - p 



a/4 + p 2 - p 



) 



,(7) 



where p = L x /L z is the aspect ratio of the stripe. The above function, which 
gives the gain in dipolar energy from the interaction of the two FM regions 
with opposite magnetization, varies between AE^-i = for p = (infinitely 
elongated sample) to AE^-i ~ — 1.9Q for p = 1 (square sample). 

The two-domain state can be favored with respect to the FM state only if 
AE dip ~ AE 0fi + A£i_i < 0. Since AE , ~ 4fi and A£i_i « -2Q, the 
conclusion is that a two-domain state with an extended domain wall can not 
appear in a one-monolayer thick stripe, whatever the aspect ratio is. For the 
triangular lattice one has AEq^ ps 5f2 and AEi,-i ~ — 2f2, so the same con- 
clusion can be drawn. This is not in contradiction with what we found in Ref. 
[8], because there we assumed an atomically sharp domain wall {Ldw = in 
the present notations). 

As stressed in the Introduction, the reason why we are now considering L^w 
1, is because exchange interaction widens the domain wall. It is well known 
[13] that in the presence of a short range ferromagnetic interaction, A, and of 
a uniaxial anisotropy, K\ (favoring the ±i directions), the additional contri- 
bution to the domain wall energy per unit length is AE ex+ani m \J AK\ > 



and the resulting domain wall width is L dw m yjA/Ki. Therefore, the total 
energy difference (per unit length) AE tota i between the two-domain and the 
FM single domain state is 



Since the gain in energy of a two-domain state must be provided by dipolar 
energy, a positive value for (A_E ,o + AEi-i) can not be healed by the addi- 
tional positive term AE ex+ani . The conclusion of this Section is that a single 
FM monolayer does not split into two domains if exchange interaction makes 
the domain wall width finite. 



3 Micromagnetic calculation of the magnetostatic energy of a rect- 
angular prism 

A few years ago, Aharoni [10] gave an analytic expression for the demagnetiz- 
ing factors of a uniform and homogeneous ferromagnetic particle in the shape 
of a rectangular prism, extending over the volume —a < x < a, —b<y<b 




AE total = A£ ,o + A£i _i + AE, 



'ex+ani- 



(8) 
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Fig. 2. Main: The dipolar contribution to the domain wall energy per unit length, 
A£^o,0) as a function of 2L Z , for 2Ld w = 5 (circles) and IL^m = 20 (diamonds). Inset: 
The same quantity, A£^o,o ; as a function of L^ w , for L z = 10 3 . All results refer to a 
square lattice. 

and —c<z<c, with the origin of the Cartesian coordinate system at the 
center of the prism (see Fig. 3, top). With respect to the notations used in 
the previous Section, we have a = a L x , b = a L y , c = a L z , and w = a L dw 
(see footnote 2 ). Here we briefly sketch Aharoni's calculation [10] for clarity's 
sake. 

When the prism is saturated along +z, surface charges ±M S are created on 
the faces z = ±c, where M s is the saturation magnetization. The potential 
due to the surface charges is (in cgs units) [13] 

s 1 1 



The density of the surface charges is (n ■ M) where M is the magnetization 
and n, the unit normal to the surface S of the ferromagnetic body, is taken 
to be positive in the outward direction. The magnetic field generated by the 
surface charge distribution is H = — VU and the magnetostatic self-energy is 



Since L^ w <C L x , the definitions o — clqL x and o — q>o(L x + L^w) are equivalent. 
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Fig. 3. The geometry of the rectangular ferromagnetic prism. The bold arrows denote 
the magnetization directions. Schematic representations of the surface magneto- 
static charges are displayed for the single-domain state (top) and for the two-domain 
state (bottom). 

given by 

E M = - X - Jm ■ H dV, (10) 

v 

where the integration is over the volume V of the ferromagnetic body. As 
a general rule [13,17], it is not advisable to evaluate first the potential U 
for a specific case, then obtain the magnetostatic field H = — V£7, and finally 
substitute H in Eq. (10). It is better to write (10) with all its integrals and then 
consider which to perform first. In other words, the order in which integrations 
are carried can considerably change the amount of algebra. The magnetostatic 
self-energy of the uniform rectangular prism depicted in Fig. 3 (top) is 

I +c b b a a 

E M = j dz Q- z J dy J d v J dx J d£ 

— c —b —6 —a —a 
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y/(x - O 2 + (y- V) 2 + (z- cY ^{x - O 2 + (y- V) 2 + + cf 

By carrying out analytically the four integrations over x, y, £, 77 (see Appendix 
A for details) one obtains 

Eti = \M 2 S VN z , (11) 



where the explicit expression for N z , the demagnetizing factor of the general 
rectangular prism of volume V = 8abc, is given in Ref. [10]. In general one has 
N x + N y + N z = 4vr; for a cubic sample (a = b = c), N x = N y = N z = ^. 

When the rectangular prism is divided into two domains with opposite mag- 
netizations directed along ±z (see Fig. 3, bottom), the contribution of surface 
charges to the magnetostatic energy is 



+c „ b b 

a 



a —w 



--c -b 

1 



-b —a —a 



y/(x - o 2 + (y - v) 2 + (z- c) 2 y/( x - o 2 + (y- v) 2 + 5 + ^ 

+c 6 6 a a 

— y ^z— J dy J dr] J dx J dt; 

1 



-6 -6 -a w 
1 



^ - O 2 + (S - v) 2 + (z- cf ^{x - o 2 + (y - v) 2 + (z + c) 



Taking into account the symmetry properties of the multiple integrals (see 
Appendix A), one obtains = 0. Thus, the difference AE M = E^ — E 1 ]] 
between the surface magnetostatic energy of the prism in the two-domain 
configuration and the magnetostatic energy of the same prism in the single- 
domain configuration is given by 

AE M = E$- E^ = -\m 2 s VN z , (12) 



i. e., it is just the opposite of the magnetostatic self-energy of the single-domain 
prism. The explicit expression of AEm, as a function of the reduced thickness 
t = b/c and of the in-plane aspect ratio p = a/c of the rectangular prism, 
turns out to be 
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^((p 2 + t 2 )f + (l+p 2 )f + (l + t 2 )f) . (13) 



Clearly, the quantity AEm (see Eq. 12) is always negative: as for the surface 
contribution to the magnetostatic self-energy, the system would prefer to di- 
vide in two domains. In the limit of vanishing thickness, the magnetostatic 
energy difference AE M approaches 0, since for t — > one has 



N z ^ -4* bit 

+ ^4 ln(2p) -1 + 2 + ^1+^ + 2 ln[l + 2p 2 yfl + p 2 } j + 0(t 2 ). (14) 

The vanishing of AEm could be expected on general grounds because the 
number of surface charges tends to as the height of the prism shrinks [18] . 

A limit to the splitting of the sample in two domains is posed by the full 
cost of the domain wall energy, which involves other energy terms: exchange 
interaction and anisotropy. The first term (exchange) is short ranged and 
favours parallel alignment of neighboring spins: 

E ex = ^ /(V • M) 2 dV>. (15) 

s v 



The anisotropy term, assumed to favor the alignment of the magnetization 
along the z axis, has the form 



F Kl 



M 2 



J {M z f dV. (16) 



The magnetostatic contribution to the Neel domain wall energy is determined 
assuming a one-dimensional model of the wall [13], i.e., the magnetization 



11 



within the wall (—w < x < w) is a function of x only. The components of the 
unit magnetization are assumed to be 



/ \ w2 t \ t \ x\/2w 2 + X 2 

m x {x) = -j- — -, m y (x) = 0, m z (x) = — — — , (17) 

W 2 + X 2 W 2 + X 2 



where w determines the wall width. At the ends of the wall, where the domains 
begin, one has m 2 (±oo) = ±1. In a thin film, the magnetostatic energy of the 
volume charges within a Neel wall which occupies the region — w < x < w can 
be approximated by [13] 



E m = AnM 2 s b J dx[m x (x)Y 



— oo 
oo 



oo oo 



-2M S 2 J dt- — I J dx J dx' cos[(x - x')t] m x (x)m x (x'), (18) 



— oo — oo 



provided that the domain wall width is much smaller than the lateral dimen- 
sions of the sample, w a < c. In the same approximation, the exchange and 
the uniaxial anisotropy contributions to the domain wall energy are, respec- 
tively 



b oo 

E ex = A j dy J dx 



dm 7 x 



dx 



En 



-b — oo 
b oo 

= Ki J dy J dx m x . 

—b — oo 



dm 2 
dx 



(19) 



One finally obtains the following expressions for the magnetostatic, the ex- 
change and the anisotropy contributions 



le 



E 

J -'m 

E ex _ 2tcA 
4bc w 

E u 7TW 



-ln(l + -) 

h w' 



(V2-1), 



lam Abe 2 u 



(20) 

(21) 
(22) 



to the total Neel wall energy = 1m + lex + 7am per unit wall area (46c). 
The actual domain wall width is determined by minimizing 77V with respect 
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to w, thus leading to the transcendental equation [13] 



2A(V2-l) K x 



2w . / b\ 
1 + — lull + - 

b v w J 



w 



= 0. 



(23) 



vrM> 2 2nM 2 s 



w + b 



We have solved Eq. (23) using the material parameters of Permalloy (i.e., 
exchange constant A = 1CT 6 erg/cm, uniaxial anisotropy constant K\ = 10 3 
erg/cm 3 , saturation magnetization M s = 800 emu). 

In Fig. 4 (top) we plot the energy density (i.e., the energy per unit wall area 
46c) versus the in-plane aspect ratio p = a/c of a rectangular prism (larger 
side 2c = 3 x 10~ 4 cm), for fixed dot thickness (2b = 3 x 10~ 6 cm). Using 
the Permalloy material parameters, the calculated energy density turns out 
to be positive for p < p c and negative for p > p c , with p c = 0.27. Therefore, 
an elongated sample (p < p c ) prefers to assume a single-domain configuration 
since the energy density cost of the Neel wall, jn, exceeds the gain in surface 
magnetostatic energy density, 5m = which is obtained splitting the 

sample into two domains. Upon increasing the in-plane aspect ratio from the 
critical value p c to p = 1 (limit of a square dot), there is an increasing gain 
in the surface magnetostatic energy, and the element prefers to assume a 
two-domain configuration. In Fig. 4 (bottom) the energy density of a square 
Permalloy prism (side 2a = 2c = 3 10~ 4 cm) is plotted as a function of the 
reduced thickness t = b/c. This means that for a square sample the two-domain 
state is always preferred, whatever is the thickness. 

These theoretical results for rectangular Permalloy prisms appear to account 
for recent MFM data [3] in a series of elliptic patterned Permalloy dots with 
thickness 2b = 30 nm. In these samples, one axis of the elliptic elements was 
kept the same (2c = 3 /im) and the other was varied, so that a series of dots 
with in-plane aspect ratio p = a/c ranging between p = 0.13 and p = 1 was 
obtained. Moreover, the interdot distance was kept sufficiently large so as to 
consider the dots as non interacting. The MFM images clearly showed that 
elongated dots with p < 0.22 are in a single-domain configuration, while dots 
with p > 0.28 are in a two-domain configuration, in fair agreement with our 
theoretical result of a critical aspect ratio p c = 0.27. For circular dots, a vortex 
configuration turned out to be preferred [3]. 

The shape of the sample is fully determined by the in-plane aspect ratio p = 
a/c and by the reduced thickness t = b/c. In Fig. 5 we plot the critical aspect 
ratio p c as a function of t, for Permalloy samples. For not so small thicknesses 
(t ~ 0.1), the sample is always in a two-domain state except for very elongated 
systems (p < 0.1). With decreasing the thickness, the value of p c increases 
and for t < 0.0015, p c > 1, which means that the system always stays in the 
single-domain state. This is due to the additional exchange+anisotropy cost 
of the domain wall, because, as it appears from Eq. (12) and Eq. (14) (see 
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Fig. 4. Top: energy per unit wall area (46c) of a rectangular ferromagnetic Permalloy 
prism versus the in-plane aspect ratio p = a/c for fixed thickness. The dotted line 
is the cost in energy density, 77V, of a Neel wall separating two opposite domains 
(see Eqs. (20-22)), the dashed line is the gain in magnetostatic energy density, 
5m = (see Eq. (13)), and the full line is the total energy density 7^ + 6m- 

The single-domain state is energetically favored when p < p c = 0.27, while the 
two-domain state is preferred for p c < p < 1. Bottom: energy per unit wall area 
(46c) of a square (p = a/c = 1) ferromagnetic Permalloy prism versus the reduced 
thickness t = b/c. In this case, the two-domain state is always preferred. 



also the next Section), the balance of dipolar energy is always in favor of the 
two-domain state for a square prism (p — 1). 
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0.001 



t=b/c 



Fig. 5. The critical in-plane aspect ratio function of the reduced thickness t, 

for Permalloy films. A value of p c larger than 1 (as it is for very small thicknesses) 
means that the system always stays in a single domain state: this is an effect of the 
exchange+anisotropy cost of the domain wall. 



4 Changing the film thickness: discrete vs continuum approaches 



In Section 2 we have shown that a one ML thick magnetic particle does not 
undergo a splitting from a single-domain to a two-domain state with pass- 
ing from an elongated to a square sample. In Section 3, within a continuum 
approach, we have shown that such splitting does exist for a film. Therefore, 
we should be led to conclude that this transition appears at some critical 
thickness. In the following we are going to discuss this problem. 

First, we observe that the key quantities to be compared are (Ei _i — Ef M -,), 
(* o - Kf) to the discrete approach, aad (£« - k), E m in the continuum 
approach. The quantities (B^l - E™) and (£« - fiJJ) are negative and 
mean the gain in dipolar energy when two domains of opposite magnetization 
are formed. The other quantities, (£"0,0 — ^oo 1 ) anc ^ E m , are positive and 
mean the dipolar cost in forming a domain wall. A square-shaped particle can 
possibly split in two domains only if the dipolar gain is larger than the dipolar 
cost. In the following we are comparing discrete and continuum expressions of 
the dipolar gain and of the dipolar cost for a square particle. 3 



3 An elongated film (p <C 1) is surely in a single domain state. A transition to a 
two-domain state occurs if and only if a square particle is in a two-domain state. 
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4-1 The dipolar gain 



According to Eq. (12), AE M = E 1 ^ - E§ = -\M 2 VN Z . The quantity N z is a 
complicated expression, but we can confine here to p — 1 (square prism) and 
t — > (vanishing thickness). In these limits, see Eq. (14), 

AE M = -]-M 2 s VN z ~ -^M s V(-4tlnt + ct), (24) 

with c ~ 2.9. In the ML limit, see Eq. (7), we had 

-£™~-1.9Q(2L 2 ). (25) 



It is worth noting that for both approaches, in the limit of a domain wall 
width much smaller than the in-plane size of the particle, the dipolar energy 
gain does not depend on the domain wall width. 

If we consider the case of domain wall width much larger than film thickness, 
Ldw ^> L y , the result (25) for the ML can be extended to finite film thickness 
simply multiplying it by (2L y ) 2 , because there are (2L y ) stripes magnetized 
in the z direction which interact with the (2L y ) stripes magnetized in the 
opposite direction, so that 

E x -i - E™ ~ -\m(2L z )(2L y ) 2 ~ -)-M 2 s V{?>M). (26) 



Therefore, the discrete approach is directly comparable with the continuum 
one: the former gives N z m 3.8t and the latter N z m — Atlnt + 2M. 



4-2 The dipolar cost 



According to Eq. (20), the dipolar cost of the Neel wall is 

1 



E m = Abe it z M z s w 



-ln(l + -J 

h w J 



(27) 



where 2w is the domain-wall width and 2b is the thickness. In the relevant 
limit w ^> b, we get 

E m » ^M 2 Vt. (28) 
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In the discrete approach for the ML, see Eq. (4) and Fig. 2, we have the 
corresponding quantity 

E 0>0 -E™ ~nd (2L z ), (29) 

with d ss 5. Therefore, for both approaches, the dipolar domain wall energy 
does not depend on w in the limit wCa (or, equivalently, L dw <C L x ). 

Extending the discrete approach to a finite thickness, we have to distinguish 
between interacting spins belonging to the same layer and to different layers. 
There are (2L y ) planes and (2L y )(2L y — l)/2 distinct pairs of planes; within 
the usual approximation L& w ^> L y , the interaction between different planes 
within the domain wall does not depend on their distance, so that 

E 0fi - E™ ~ Q(2L Z ) 
where we have retained only the leading term. 

In the above expressions, d ~ 5, while d\ represents the interaction energy 
between different layers within the domain wall. We have numerically checked 
that d\ > 5, therefore the dipolar energy cost (di/2)M 2 Vt always dominates 
on the dipolar gain (— 1.9M 2 Vt). On the other hand, in the continuum picture 
we have a dipolar gain AE M = — ^M 2 V(— 4t lnt + 2.9t) which dominates on 
the dipolar cost E m = (ir 2 /A)M 2 Vt for any thickness t. 

To sum up, the continuum approach suggests that a square particle is splitted 
in two domains for any thickness, while the discrete approach supports the 
opposite conclusion. However, extrapolation of the discrete results from the 
ML to finite thickness is valid only within the limit of thickness much smaller 
than the domain wall width. We are therefore led to conclude that a critical 
thickness L* ss L^w should exist, such that the transition between single- 
domain to two-domain state appears at thicknesses L y > L* only. 

5 Conclusions 

In Ref. [8] we had considered a single monolayer with planar spins interact- 
ing only via dipolar interaction. For a triangular lattice, it was possible to 
observe a transition from a single-domain to a two-domain to a macrovor- 
tex state, when passing from elongated to square (or circular) samples. This 
phenomenology did not apply to a square lattice, because in that case spins 
do not like ferromagnetic configurations. The domain wall in the two-domain 
state was atomically sharp. 



d (2L y ) + d 1 



(2L y )(2L y -l) 



^M*Vt, 



(30) 
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In this paper we have studied in detail the transition from a single-domain to a 
two-domain configuration for a film of general thickness and in the presence of 
exchange interaction and single ion anisotropy. The presence of the exchange 
interaction has a twofold effect: it widens the domain wall and, together with 
the anisotropy, it increases the energetic cost of the domain wall. The first 
effect (a larger domain wall) is enough to suppress the transition from single 
domain to two-domain state in a monolayer and in ultrathin films. These 
results have been obtained within a discrete approach. 

For thicker films a continuum, micromagnetic approach is more appropriate. 
According to it, the energetic balance of the sole dipolar interaction is such 
that the above transition always occurs on passing from elongated to square 
samples. When exchange and anisotropy are correctly taken into account in 
the full energy balance, it is straightforward to realize that the transition may 
only occur for systems where the anisotropy K\ is vanishing small, because the 
exchange A is usually so large that the additional cost in domain wall energy 
density \J AK\ can not be compensated by the gain in dipolar energy density, 
of order M^b. A material satisfying the condition \fA~K[ <C M^b is Permalloy. 
For rectangular prisms made of this material, we have therefore determined 
the critical in-plane aspect ratio function of the prism thickness. Fair 

agreement with existing experimental data on Permalloy particles [3] was ob- 
tained. 

A final comment concerns the discrete and the continuum approaches. The 
former has been applied to ultrathin films and has shown that increasing the 
domain- wall width suppresses the transition. In contrast, the transition is ob- 
tained using the latter approach, valid for thicker films. Therefore, these results 
suggest the existence of some critical thickness d* above which the transition 
appears. The discrete formalism suggests that this thickness d* should be of 
order of the domain wall width w, but the continuum formalism does not seem 
to indicate that something "critical" should occur when b ~ w. This is the 
sole issue we have not been able to clarify in full. 
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A Analytical calculation of the multiple integrals 



The magnetostatic self-energy of the distributions of surface charges depicted 
in Fig. 3 can be expressed in terms of the multiple integrals 



+a a 



I + (w,a,b,k) = J dy J drj J dx J d£ 



-b —a w 

+a —w 



I (w, a, b, k) = J dy j drj J dx J dt; 



—b -b —a 



I y/(x - 2 + (y - V) 2 + k 2 



(A.l) 



By symmetry, it immediately turns out that I + (w,a,b,k) = I~(w,a,b,k) = 
I(w,a,b,k) and that I(w,a,b, —k) = I(w,a,b,k). The explicit expression of 
I(w, a, b, k) is 



I(a, 6, k) = l(k 2 )l - \{Aa 2 + Ab 2 + k 2 )l + \{Aa 2 + k 2 )l + \{Ab 2 + k 2 )l 



((w + a) 2 + k 2 )i - ((w + a) 2 + Ab 2 + k 2 )^ 
- ((w - a) 2 + k 2 )l + ((w - a) 2 + Ab 2 + k 2 )l 



k 2 VAa 2 + Ab 2 + k 2 - k 2 VAa 2 + k 2 - k 2 VW+¥ 



+ k 2 



^J(w + a) 2 + k 2 - ^J(w + a) 2 + Ab 2 + k 2 



- sj{w - a) 2 + k 2 + ^(w - a) 2 + Ab 2 + k 2 

o,, ( 4a6 \ 

— Sabk arctan 



+ Abk 



(w + a) arctan 
— (w — a) arctan 
+ bk 2 In 



k\/Aa 2 + Ab 2 + k 2 J 

2(w + a)b 



M^(w + a) 2 + Ab 2 + k 2 , 
2(w-a)b V 
kJ(w-a) 2 + Ab 2 + k 2 ). 



, , , VW+W + 2b" 



VW+k 2 -2b / 



- bk z In 



o , , VAa 2 + Ab 2 + k 2 + 26 s 



VAa 2 + Ab 2 + k 2 - 2b, 



+ Aa z b\n 



, ; , , V4a 2 + Ab 2 + k 2 + 2b\ 



^Aa 2 + Ab 2 + k 2 -2b) 
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+ 2ak 2 In (^/Aa 2 + k 2 + 2aj 
+ k 2 -(w + a) In (^J(w + a) 2 + k 2 - (w + a)^J 
+ (w — a) In ( \J ' (w — a) 2 + k 2 — (w — a) 



- 2ak 2 In \^/Aa 2 + 4b 2 + k 2 + 2a j 



+ k 2 



w + a) In J(w + a) 2 + Ab 2 + k 2 - (w + a) 



(w — a) In (w — a) 2 + Ab 2 + k 2 — (w — a) 
8ab 2 In (^/Aa 2 + Ab 2 + k 2 - 2a^j 

w + a) In I y/{w + a) 2 + Ab 2 + k 2 - {w + a) I 
(w - a) In (^J(w - a) 2 + Ab 2 + k 2 - (w - a)j 



+ b 



(k 2 -{w + a) 2 ) In 



- (k 2 -(w- a) 2 ) In 



(w + a) 2 + Ab 2 + k 2 + 2b s 
(w + a) 2 + Ab 2 + k 2 -2b, 
[w - a) 2 + Ab 2 + k 2 + 2b\ 
'w - a) 2 + Ab 2 + k 2 - 2b) 



— k 2 (w + a) In 
+ (w — a) In 
+ k 



w + a) 2 + k 2 + (w + a) 
y/(w + a) 2 + k 2 -(w + a) 
\J{w — a) 2 + k 2 — (w — a)' 
y 1 (w — a) 2 + k 2 + (w — a) , 



(w + a) In 
+ (w — a) In 



> + a) 2 + Ab 2 + k 2 + (w + a) 



[w + a) 2 + Ab 2 + k 2 -(w + a 
'w — a) 2 + Ab 2 + k 2 — (w — aY 



^J(w - a) 2 + Ab 2 + k 2 + (w - a) 



(A.2) 



Taking the limit w — > in the terms in square parentheses on the r.h.s., one 
obtains 



1(0, a, b, k) = \ ( k 2 ) § - \ (4a 2 + Ab 2 + k 2 ) § 
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AV4a 2 + 4b 2 + k 2 



2 + k 2 + 2b 



+ l(4a 2 + A; 2 )l + i(46 2 + A; 2 )f 

+ fc V4a 2 + 4fr 2 + P 

- fc 2 V4a 2 + A; 2 - A;V46 2 + P 

/ 4a6 

- 8abk arctan 
+ bk 2 In 

- bk 2 In j 
+ 4a 2 Mn 
+ ak 2 In I 



ak 2 In 



\VWTk 2 -2bJ 
V4a 2 + 46 2 + A; 2 + 2&\ 
v / 4a 2 + 46 2 T~P-26/ 
V4a 2 + 46 2 + A; 2 + 26 s 
V \/4a 2 + 46 2 TI 2 -26 / 
v^a 2 + A: 2 + 2a \ 
V4a 2 + A; 2 ^2a ) 
\/4a 2 + 46 2 + A; 2 + 2a\ 



V4a 2 + 46 2 + k 2 - 2a 



. u2 , fV4:a 2 + 4b 2 + k 2 + 2a\ . . „, 

+ 4a& 2 In (A. 3) 

Vv / 4a 2 + 46 2 + P-2a/ V ; 



Dropping the a, 6 arguments for brevity's sake, and taking into account the 
symmetry properties of the multiple integrals, the magnetostatic self-energy 
for the one-domain configuration (see Fig. 3, top) is expressed as 



= ^ I M 2 



4tt 2 



r(o, o) - r(o, -2c) + i + {o, o) - 7+(o, -2c) 



J-(0, 2c) + /-(0, 0) - J+(0, 2c) + J+(0, 0) 



1B_ 

An 



M 2 2 



7(0,0) -7(0, 2c) 



(A.4) 



Upon substituting the explicit expression of the multiple integral (A. 3) in the 
equation (A.4), it is immediate to recover Aharoni's formula for the demagne- 
tizing factor of the uniform and homogeneous ferromagnetic rectangular prism 
(Eq. 1 of Aharoni's paper) [10]. 

For the two-domain configuration (see Fig. 3, bottom), the contribution of the 
surface charges to the magnetostatic energy is found to be 



E 



IT 

M 



4tt 2 3 



r(w, 0) - r(w, -2c) - I + (w, 0) + I + (w, -2c) 



r(w, 2c) + I~(w, 0) + I + (w, 2c) - I + (w, 0) 



(A.5) 



so that, taking into account the symmetry properties of the integrals, E 



M 
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is found to vanish. In conclusion, the magnetostatic energy difference AEm 
turns out to be just the opposite of the magnetostatic self-energy of the single- 
domain prism. Its explicit expression is given in the main text, Eq. (13), in 
terms of the in-plane aspect ratio, p = a/c, and of the reduced thickness, 
t = b/c, of the rectangular prism. 
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